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We present a systematical study on the prompt J/ip production in association with a cc pair 
via the process, 77 —> H (cc) + c + c, within the framework of non-relativistic QCD at the future 
high-energy e + e _ collider - International Linear Collider (ILC), including both direct and feed-down 
contributions. For direct J/ip production, the states with color-octet channels, especially the 3 Pj^ 
and Sq j ones, provide dominant contribution to the production cross-section, which are about 
fifty-two times over that of the color-singlet one. This is clearly shown by the transverse momentum 
(pt) and rapidity distributions. The feed-down contribution from if' and Xcj (J = 0 , 1 , 2 ) is sizable, 
which is about 20% to the total prompt cross-section. Besides the yields, we also calculate the J/ip 
polarization parameter A. In small pt region, the polarization of the prompt J/ip is longitudinal 
due to large contributions through the 3 Pj^ channel, and becomes transverse in high pt region due 
to 3 S® channel. Thus the J/ip production via photon-photon collisions at the ILC shall provide a 
useful platform for testing the color-octet mechanism. 

PACS numbers: 13.66.Be, 12.39. Jh, 14.40.Pq, 12.38.Bx 


I. INTRODUCTION 

With its heavy mass m the heavy quark (Q) will move 
with a small velocity vq inside a heavy quarkonium. This 
results in a hierarchy of energy scales, m » mvQ » 
mvQ, and the dynamics at different energy scales is dif¬ 
ferent. The non-relativistic QCD (NRQCD) provides a 
systematical way to separate the effects from the dynam¬ 
ics at different energy scales [1], Since its invention, 
it has been widely applied to deal with heavy quarko¬ 
nium physics. It provides a solution to the if' {J/ip) 
surplus puzzle [2-4] . It gives a consistent explanation for 
both the J/ip polarization and the r) c production data 
at hadron colliders [5-7] by using the newly fitted color- 
octet long distance matrix elements (LDMEs) [8, 9]. It 
gives an explanation for the x c j-production data mea¬ 
sured by the hadronic colliders [10]. A comprehensive 
review of NRQCD theory and applications can be found 
in Refs. [11, 12], 

The NRQCD theory still faces many challenges. As an 
important example, the Belle and the LHCb collabora¬ 
tions [13, 14] have measured the J/ip production associ¬ 
ated with a cc pair, which however show large discrepan¬ 
cies from the NRQCD predictions [15-22]. For instance, 
the predicted angular distribution of J/ip is different from 


the Belle data, and the predicted J/ip production cross- 
section is lower than the LHCb data by about one order 
of magnitude. This discrepancy may indicate a breaking 
of NRQCD factorization in small J/ip momentum region 
at the Belle. Thus it is helpful to find more experimen¬ 
tal platforms to test NRQCD, especially those with high 
collision energies. 

A detailed prediction for the J/ip production associ¬ 
ated with a cc pair, e + e _ —> J/ip+cc , at the super-Z fac¬ 
tory has been done in Ref. [23] . Due to the Z°-boson res¬ 
onance effect, sizable J/ip events can be achieved there. 
The proposed International Linear Collider (ILC) [24] is 
another useful e + e _ collider, which has been designed 
to run at a high collision energy from several hundred 
GeV to TeV together with a high luminosity about 
£ ~ 10 34 ” 36 cnU 2 s _1 . An analysis of J/ip production in 
association with a cc pair via photon-photon collision at 
the ILC has been done by Ref. [25]. It shows sizable J/ip- 
events, about 4.7 x 10 4 per year, can be generated under 
the condition of \/S = 500 GeV and £ ~ 10 34 cm _2 s _1 . 
As a comparison, the J/ip photon-photon production in 
association with light hadrons is very small [26] 1 . 

It is noted that in Ref. [25], the J/ip is generated via the 
direct production channel, 77 —¥ J/ip + cc., in which only 
the color-singlet contribution has been taken into consid¬ 
eration. According to our NRQCD practices, the contri- 
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1 It only contributes 0.4% and 0.02% to the J/tp inclusive produc¬ 
tion cross-section with y/S = 500 GeV and 1 TeV, respectively. 
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buttons from the color-octet charmonium states may be 
sizable or even dominant. Furthermore, the feed-down 
contributions from the ip' and XcJ may also be sizable, 
i.e. one can first generate the ip' and XcJ mesons via the 
channel, e + e _ —> ip'/XcJ + cc, and then they may decay 
to J/ip via the channels ip' -A J/ip + X , y c j -A J/ip + 7 
and ip' -A XcJ —> J/ip + 7 - Thus, as a sound a NRQCD 
prediction, we shall in this paper study the prompt J/ip 
production with a cc pair, including both the direct and 
the feed-down contributions from ip' and XcJ mesons. 

For the direct J/ip and ip' production, four interme¬ 
diate states shall be considered, i.e. the color-singlet 
one ( 3 S[ 1] ) an d the three color-octet ones ( 1 5'g 8 ^, 3 sj 8 ^ 
and 3 P^), which keep the NRQCD expansion up to the 
order of Vq. For the direct XcJ production, up to the 

leading order in vq, only the two dominant states ( 3 P 
and 3 S'[ 8 ^) shall be considered. The feeddown contribu¬ 
tions from ip' to XcJ shall also be considered in our work, 
which, as shall be shown, give ~ 2 % contribution to the 
prompt J/ip production. 

The remaining parts of the paper are organized as fol¬ 
lows. In Sec.II, we describe the calculation technology 
for treating the J/ip production in association with a cc- 
pair. The input parameters are also introduced there. In 
Sec.Ill, we give the predictions for the yields and polar¬ 
izations of J/ip at the ILC. The last section is reserved 
for a summary. 


II. CALCULATION TECHNOLOGY AND 
INPUT PARAMETERS 

At the ILC, the charmonium H can be produced via 
the process e + + e~ -A e + + e~ + H(cc) + c + c. For 
the photon-photon production, 77 — > H(cc ) + c + c, 
the photons may interact either directly (direct photon 
production) or indirectly via their hadronic components 
(resolved photon production) with the quarks or anti¬ 
quarks. In the paper we shall focus our attention on 
leading contribution from the direction photon produc¬ 
tion. Within the NRQCD framework, the differential 
cross-section for the charmonium H can be factorized as 

da H = / dx 1 dx 2 f 1 {xi)f 1 (x 2 ) 

n ^ 

x dd (77 —> cc[n] + c + c)(O h (n)). ( 1 ) 

a stands for short-distance cross-section, representing the 
production of an intermediate perturbative state ( cc)[n] 
with quantum number n. (O h (n)) is non-perturbative 
but universal LDMEs, which cannot be evaluated per- 
turbatively, but can be obtained through the fit of the 
experiment. f-y[x) is the photon density function with x 
being the momentum fraction of the photon to the initial 
electron or positron. 

For the hard subprocess 77 -A cc[n] + c + c, there are 
totally 20 diagrams for [n] = 3 5'[ 1 ^ and [n] = 3 pj^, 32 dia¬ 



FIG. 1: Typical Feynman Diagrams for the hard subprocess 
7(Pib(P2) -A cc[n](p 3 ) + c(p4,) + c(p t ), in which (o, d ) are for 
color-singlet states, (a, b, c, d) are for cc[ 3 S® ] and (a, d, e) are 
for cc[ 1 5q 8 ^] and cc[ 3 Pj 8 '], respectively. 


grams for [n] = 3 S'| 8 ^, and 28 diagrams for [n] = 1 S'q 8 ^ and 
[n] = 3 P^. Typical Feynman diagrams are presented 
in Fig.(l). Those Feynman diagrams can be generated 
and the corresponding short-distance cross-sections can 
be calculated by using the well-established FDC package. 
The FDC is a general-purpose program for Feynman di¬ 
agram calculation which realizes an automatic deduction 
from physical models to final numerical results for lower- 
order processes [27]. 

At the high-energy e + e _ collider such as ILC, the 
laser backscattering (LBS) from the incident electron 
and positron beams leads to high luminosity photon 
beams [28]. Those LBS photons are hard enough and 
carry a large fraction of the energy of the incident lepton 
beams. Moreover, as pointed out by Refs. [26, 29, 30], 
using LBS photons the 77 collisions can approximately 
achieve the same high luminosity as that of the e + e“ 
beams. Thus, we adopt the LBS density function to do 
our calculation [28] 


fy(x) 


1 

N 


1 — x + 


1 


1 — X 


4r(l — r) , 


( 2 ) 


where r = x/[x m {\ — x)\ and the normalization factor 




logx 



( 3 ) 


Here y = 1 + x m , x m ~ 4.83 [31] and the LBS photon 
energy is restricted by 0 < x < x m /{l + x m ) = 0.83. 

The unpolarized differential cross-section of prompt 
J/ip is obtained by adding the unpolarized cross-sections 
of the various direct-production processes multiplied by 
appropriate branching fractions, i.e. 

^.prompt J/d 

= da J ^ BrpXcJ ~► J/ip + 7 ) dcr XcJ + 

J 

y: Brpip' -A XcJ + 7 )Br(xcj —t J/ip + 7 ) da ^ 

J 

+Br(iP' -A J/ip + X) da*'. 


( 4 ) 





3 


As for the branching ratios, we take [32] 


-A . 

J/ip + X) 

= 60.9%, 

Br{ip' — 

i XcO + 7) 

= 9.99%, 

Br(ip' - 

► Xci + 7) 

= 9.55%, 

Br(i/)' - 

I Xc2 + 7) 

= 9.11%, 

Br(xc0 -a 

J/i/j + j) 

= 1.27%, 

#KXcl -a 

J/i/j + j) 

= 33.9%, 

Br(xc2 -A 

J/ip + x) 

= 19.2%. 

As for the LDMEs, ’ 

we take [8- 

-10]: 


2 . Regarding the feed-down contributions, the resul¬ 
tant J/t/’ transverse momentum can be achieved by a 
shift, = "If/* Pt 1 , where to# stands for charmo- 

nium mass. The charmonium masses are [32]: mj= 
3.097 GeV, m r = 3.686 GeV, m Xc0 = 3.415 GeV, 
m Xcl = 3.511 GeV, and m Xc2 = 3.556 GeV. 


VS 


a 


J/i> 

3S] 1 ] 


250 GeV 9.28 x 10“ 4 
500 GeV 4.03 x 10" 4 
1000 GeV 1.58 x 10“ 4 


1 S [8] 

_ 


1.40 x 10" 3 
1.09 x 10" 3 
5.66 x 10“ 4 


<7 


J/ip 

3S{ 81 


3.46 x 10“ 5 
1.62 x 10“ 5 


6.87 x 10“ 6 


a 


3 P] 81 


1.94 x 10" 2 
1.46 x 10" 2 
7.60 x 10“ 3 


{ oxc J( 3 P m )} 

(o^'es®)) 

(0^{ 3 s/ ] )) 

ml 

(o^'Csf/) 

(O j / ,/ '( 1 1 S'' 81 )} 

(O j ^(/S 1 i ] )) 

(O j ^{ 3 p1 8] )) 


(2J+ 1) x 2.15 x 10" 3 GeV 3 , 
(2J + l)x A|i?' 1P (0)| 2 , 

1.0 x 10~ 2 GeV 3 , 

0.401 x 10" 2 GeV 3 , 

0.682 x 10“ 2 GeV 3 , 

0.645 GeV 3 , 

0.780 x 10“ 2 GeV 3 , 

1.057 x 10" 2 GeV 3 , 

1.934 x 10“ 2 GeV 3 , 


where J = (0,1, 2), the squared first derivative of the ra¬ 
dial wave-function at the origin |i?' 1P (0) | 2 = 0.075 GeV 5 
and the squared radial wave-function at the origin 
|i?2s(0)| 2 = 0.529 GeV 3 [33]. We have taken the cen¬ 
tral value of (0 J /^( 3 sf^)} to be the one determined 
from a recently determined (O r,c ( 1 s/ 1 /) [8], which is 
different from the one extracted from measurements of 
T(J/ip —> l + l~). In Ref.[8], the color-singlet LDME for rf c 
has been determined via a direct fitting of the new LHCb 
data on the rj c production [7], which gives, ( 0 Vc ( 1 S'g 1 ^)} = 
0.215 ± 0.135 GeV 3 . This value is comparable with the 
ones obtained by other groups, e.g., 0.39 GeV 3 [33] and 
0.437^0^5 GeV 3 [34], By further using the heavy-quark 

spin symmetry, (0 Vc ( 1 5'q”^)} = ^(O j ^( 3 s//), we ob¬ 
tain (CV/^sj 11 )) = 0.645 ± 0.405 GeV 3 . 


III. NUMERICAL RESULTS AND 
DISCUSSIONS 


TABLE I: The integrated cross-section (in unit: nb) for the 
direct J/i/) photon-photon production in association with a cc 
pair at the ILC. \y\ < 4.5. 


Vs 

^direct 

^feed-down from i// 

^feed-down from XcJ 

250 GeV 

2.18 x 10 -2 

5.79 x 10" 3 

6.29 x 10 -4 

500 GeV 

1.61 x 10“ 2 

4.18 x 10” 3 

4.16 x 10" 4 

1000 GeV 8.33 x 10" 3 

2.14 x 10 -3 

2.07 x 10" 4 


TABLE II: The integrated cross-section (in unit: nb) for the 
prompt J/ij) photon-photon production in association with 
a cc pair at the ILC. As for the direct J/ij) production, the 
contributions through the four channels have been summed 
up. \y\ < 4.5. 

We present the integrated cross-section for the J/ij; 
photon-photon production in association with a cc pair 
in Tables I and II. The integrated cross-section decreases 
with the increment of the e + e _ collision energy y/S. The 
photon-photon production of J/ip is dominated by direct 
production, while the feed-down contribution is sizable. 
More explicitly, when \/S = 1 TeV, the direct channel 
provide ~ 78% contribution to the prompt J/ip produc¬ 
tion, that of the feed-down from if'/ is ~ 20%, and that 
of feed-down from \cJ is ~ 2%. Furthermore, for the di¬ 
rect production, the color-octet contributions are large; 
two channels ( 3 P}/ and 1 S'o 8 ^) provide over 95% contri¬ 
bution to the direct J/ij] cross-section. Thus the future 
new measurements on J/ip photoproduction at the ILC 
can provide a useful platform for testing the NRQCD 
color-octet mechanism. 

As an important point, we provide an explanation of 
the reason why 1 S'q 8 ^ and 3 pf^ channels provide such 
large contributions in comparison with the other two. 
By analyzing the topologies of Feynman diagrams, it is 
noted that the largest contributions of 1 S'q 8 ^ and 3 pj 8 ^ 
mainly come from Fig.(l.(e)), which is absent in the 


To do the numerical calculation, we take m c = 1.5 GeV 
and a = 1/137. The one-loop a s running is used, and 
for each charmonium H , we set its renormalization and 
factorization scale as pa = p,F = M t = i/4m[? + ( p/) 2 


2 Taking //// = M/,/2, the prompt cross-section shall be increased 

by 70%, and taking fiR = 2 Mt, the prompt cross-section shall be 

decreased by 30%. Such a large scale error could be suppressed 

by a high-order calculation and/or a proper scale setting [35]. 
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FIG. 2: The LBS photon density function. 


other two channels due to the C-parity and color con¬ 
servation. The squared invariant mass of the internal 
gluon of Fig.(l.(e)) is, k 2 = (pi - p 3 ) 2 = 4m;? - 2pi • p 3 
with pi-p 3 = x^f(E H -p 3z ), where E H (= e " + 2 e " M t ) is 
the charmonium H energy, p 3z (= 6 ~ e M t ) is the pro¬ 
jection of its 3-momentum to the flying direction of the 
initial photon which attaches to iT, and y is the rapidity 
of H. H stands for 1 Sq 8 ^ or 3 P^, and x is the momentum 
fraction of this photon to the electron or positron. Then 
the fc 2 -expression can be rewritten as 

k 2 = 4m^ — xVSM t e~ v 

= 4m\ - 2 [E h + E g ) x M t e~ v 
= -4m 2 c e- 2y - (1 + e~ 2y )(pf) 2 - 2E g M t e~ y , (5) 

where E g is the energy of the internal gluon. In small pt 
region, the magnitude of k 2 reduces exponentially with 
the increment of y and could be very small, leading to 
large contributions to the production cross-section. 

In addition to the topology, the LBS photon density 
function provides another important factor for large pro¬ 
duction cross-sections. A bigger energy Eh corresponds 
to a bigger rapidity \y\ 1 and vice versa. On the other 
hand, a larger photon momentum fraction x indicates 
more H events with larger energy can be achieved. As 
shown by Fig. (2), the LBS photon density function be¬ 
haves moderately in small and intermediate x regions and 
increases fast when approaching to its allowable largest 
x. Thus the contributions of Fig.(l.(e)) will be further 
enhanced by LBS photons. 

We present the pt-distributions of the direct and 
prompt J/ip productions associated with a cc pair at the 
ILC with \/S = 1 TeV in Figs. (3, 4). Fig. (3) shows that 
the color-octet and color-singlet channels have different 
p t behaviors. Fig. (4) shows the relative importance of 
the feed-down contributions from ip' and XcJ to the di¬ 
rect one. The fact that the 3 pj 8 ^ channel, together with 



FIG. 3: The pt-distributions for the direct J/ip production as¬ 
sociated with a cc pair at the ILC. The dash-dot, the dashed, 
the dotted and the solid lines are for 3 s[^, 1 S/ 8 / 3 S'® and 

q [Q] 

Pj , respectively. 



FIG. 4: The pt-distributions for the prompt J/ij) production 
associated with a cc pair at the ILC. The dashed line is for 
direct J/ip production, in which the contributions through 
the four channels have been summed up. The dash-dot line 
is for the feed-down contribution from ip' and XcJ decay. The 
solid line is for the prompt production. 


1 S'g 8 ^ channel, provide the dominant contributions to J/ip 
production cross-section can be explained by their large 
contributions in small p t region. More explicitly, the 1 S'® 
Pt distribution is over the color-singlet one for p t < 20 
GeV, while the 3 P^ distribution dominates the color- 
singlet one for pt < 45 GeV. Different from other chan¬ 
nels which scale at least as 1/p®, the 3 S'[ 8 ^ channel scales 
as 1/p/ due to the one gluon fragmentation mechanism. 
Thus even though its differential cross-section is small in 
low p t region, it becomes important in high p t region, 
i.e., p t > 50 GeV. 
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For our present considered production channel 77 —> 
J/ip + cc, there is no ln(p t /m c )-type large log-terms 
at the tree-level. This is because all the final par¬ 
ticles are massive and the intermediate gluon should 
also be hard enough either to generate a cc-quark pair 
[Figs. (la,lb, lc)] or to pull the c and c in different z di¬ 
rection together to form a J/ip [Fig. (Id)] or to put the 
free c (or c) changing its 2 -direction to run along with 
the same z-direction of J/ip [Fig.(le)]. More explicitly, 
for the dominant 1 S'q 8 ^ and 3 P^ channels, Eq.(5) in¬ 
dicates that the momentum of the intermediate gluon 
can be small but cannot equal to zero, thus the J/ip pt~ 
distributions are under well control. This can be seen 
more clearly from Figs. (3, 4) that the maximum differen¬ 
tial cross-sections appear at about pt ~ m c , and it drops 
down for even smaller p t values; at pt = 0 , its cross- 
section is suppressed from the phase-space due to all the 
final particles are massive. For a higher-order calculation, 
there may emerge ln(pj/m^)-type large log-terms due to 
extra soft parton (gluon or light quarks) coming into con¬ 
tribution, requiring special treatment, which however is 
out of the range of the present paper 3 . 



y 


FIG. 5: The ^-distributions for direct J/ip production associ¬ 
ated with a cc pair at the ILC. The dash-dot, the dashed, the 
dotted and the solid lines are for 3 sf ] , 1 S',], 81 , 3 s[ 8] and 3 Pj 8] , 
respectively. 


We present the y-distributions for the direct and 
prompt J/ip productions associated with a cc pair at the 
ILC with = 1 TeV in Figs. (5, 6 ). Fig. (5) shows 


3 For the present considered photon-photon collision processes 
with massive final particles, in different to the case of hadronic 
productions [36, 37], there is no large log-contributions from 
initial-state gluon showers. The gluon radiation is also sup¬ 
pressed from the final heavy quark lines [38, 39]. Thus the higher 
fixed-order pQCD prediction may still be trustable in small pt 
region, at least for pt > 2 m c . This could be similar to a NLO 
calculation of e+e — -4 J/ip + cc at the B factories [22]. 



y 

FIG. 6: The y-distributions for the prompt J/ip production 
associated with a cc pair at the ILC. The dashed line is for 
direct J/ip production, in which the contributions through 
the four channels have been summed up. The dash-dot line 
is for the feed-down contribution from ip 1 and y c j decay. The 
solid line is for the prompt production. 


that the rapidity distributions of the four channels can 
be divided into two groups, whose behavior is either con¬ 
cave or convex. Those two types of rapidity distributions 
can be adopted to distinguish the color-octet states from 
the color-singlet one. The concave behavior of 1 5'q 8 ^ and 
3 pj 8 ^ channels is consistent with the explanation (5) for 

the large cross-section of 1 S'q 8 ^ and 3 Pj 8 ^: the largest k 2 
is achieved at y = 0 which leads to the lowest point of 
the concave, larger rapidity leads to smaller k 2 and larger 
differential cross sections. 



y 


FIG. 7: The y-distributions under various pt cuts for the 
prompt J/ip production associated with a cc pair at the ILC. 
The solid, the dotted, the dashed and the dash-dot lines are 
for pt > 0 GeV, 2 GeV, 4 GeV, and 6 GeV, respectively. 
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contribution from ip'. In fact, because the \ c j feed-down 
contribution to the integrated cross-section is only about 
2%, its contribution to A is negligible. It is noted that 
the 3 S'[ 1 ^ polarization and the prompt J/ip polarization 
behave quite differently from each other in whole pt 
region. In small p t region, the prompt polarization 
is dominated by 3 pj 8 ^ and is longitudinal. In high p t 
region, it becomes transverse due to the one gluon 
fragmentation mechanism of 3 S^. By contrast, the 
polarization of the color-singlet channel is slightly trans¬ 
verse in small pt region, while it is almost unpolarized 
when p t > 20 GeV. Thus, besides the aforementioned p t - 
and y- distributions, the polarization parameter A can 
be another useful tool to test the color-octet mechanism. 


FIG. 8: J/ip polarization parameter A as a function of pt- 
The label ‘p rom pt’ represents the sum of the direct and ip' 

feed-down contributions. SUMMARY 


Considering the detectors’ abilities and in order to offer 
experimental references, we calculate the rapidity distri¬ 
butions under various pt cuts. The prompt J/ip rapidity 
distributions for p t > 2 GeV, pt > 4 GeV and p t > 6 
GeV are presented in Fig. (7). With the increment of the 
Pt cut, up to ~ 10 GeV, the rapidity shapes are still 
dominated by 1 S'q 8 ^ and 3 channels, following their 
concave behavior. The degree of the concavity decreases 
with with the increment of the pt cut, which is due to the 
fact that the term [—(1 + e~ 2 y )(pt) 2 ] in Eq.(5) becomes 
more and more important for a larger p t value, providing 
a suppression for the cross-sections with larger rapidi¬ 
ties. By taking a much higher p t cut, e.g. p t > 40 GeV, 
the other two direct channels, 3 S'[ 1 ^ and 3 sf\ becomes 
important, e.g. for small rapidity region \y\ ^ 1.6 their 
differential cross-sections are larger than those of 1 5 q 8 ^ 
and 3 Pf^ channels, and the prompt J/ip rapidity curve 
becomes a flat line for \y\ ^ 3. 

Reviewing the long standing J/ip polarization puz¬ 
zle which challenges the NRQCD theory [40], besides 
the yields, we make a discussion on the polarization of 
prompt J/ip of the process. The polarization observable 
A for J/ip is defined as [41] 


A = 


- da^ 4 ’ 
d<T\\’ + dcroo 0 


( 6 ) 


where dcr J JJ/,(S z ,S' z = 0,±1) is the spin density ma¬ 
trix, which can be calculated by FDC. The feed-down 
contributions from ip' and \cJ are much more involved, 
detailed procedures for calculating the parameter A for 
the prompt J/ip can be found in Refs. [41, 42]. 

We present the polarization parameter A as a function 
of pt in Fig. (8). As shown in Table II, the direct 
J/ip dominates the production cross-section, so in our 
analysis, we shall only consider the dominant feed-down 


In the paper we have studied the photon-photon pro¬ 
duction of prompt J/ip in association with a cc pair at 
the future collider ILC within the framework of NRQCD. 
The color-octet channels, especially 3 P}f > and 1 5 , q 8 ^, pro¬ 
vide dominate contributions to the production in small 
and medium p t region. 

At the ILC with the e + e _ collision energy \Z~S = 1 
TeV, the color-singlet cross-section cr 3 „[i] is only ~ 1.5% 

‘-’l 

of the NRQCD prompt prediction ctnrqcd 4 , which in¬ 
cludes both direct and feed-down contributions. The to¬ 
tal feed-down cross-section is sizable, providing ~ 22% of 
o’nrqcd, thus those feed-down channels should be taken 
into consideration as a sound prediction. Moreover, the 
predicted J/ip pt- and y- distributions, as well as the 
J/ip polarization, given by the color-singlet mechanism 
and NRQCD are quite different. 

If taking the ILC luminosity as C ~ 10 34 cm _2 s _1 , siz¬ 
able J/ip events can be generated in one operation year, 
i.e. about 2.7 x 10 6 , 2.0 x 10 6 or 1.1 x 10 6 events can be 
generated for \fS = 250 GeV, 500 GeV or 1 TeV, respec¬ 
tively. For the case of \fS = 1 TeV, the generated events 
shall be changed down to 5.0 x 10 5 , 1.5 x 10 5 , and 5.7x 10 4 
for p t > 2 GeV, p t > 4 GeV and p t > 6 GeV, respectively. 
All those are sizable quantities, the J/ip photon-photon 
production channel shall thus provide a useful platform 
for testing the NRQCD color-octet mechanism. 

One thing need to mention is that in the paper, we do 
not calculate the next-to-leading order QCD correction 
to the channels, which might be significant. However, 
due to the a s -power suppression at higher orders, we 
can expect that the contributions from the color-octet 


4 Throughout the paper, we have adopted a conservative rapidity 

range \y\ < 4.5 to do the calculation. If taking a smaller rapidity 

region as \y\ < 2.0, the percentage of cr [i] to ctnrqcd shall be 

s i 

raised up to ~ 10%. 
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channels are at least important. At the ILC, the J/ip 
can also be produced in associated with an open 56-pair, 
which is out of the range of the present paper. 
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